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Abstract. A weakly dispersive sub-range (WDR) of kinetic 
Alfven turbulence is distinguished and investigated for the 
first time in the context of MHD/kinetic turbulence transi- 
tion. We found perpendicular wavenumber spectra oc fcj 3 
and oc fc^ 4 formed in WDR by strong and weak turbulence 
of kinetic Alfven waves (KAWs), respectively. These steep 
WDR spectra connect shallower spectra in the MHD and 
strongly dispersive KAW sub-ranges, which results in a spe- 
cific double-kink (2-k) pattern often seen in observed turbu- 
lent spectra. The first kink occurs where MHD turbulence 
transforms into weakly dispersive KAW turbulence; the sec- 
ond one is between weakly and strongly dispersive KAW 
sub-ranges. Our analysis suggests that the partial turbulence 
dissipation due to amplitude-dependent super-adiabatic ion 
heating may occur in the vicinity of the first spectral kink. A 
threshold-like nature of this process results in a conditional 
selective dissipation affecting only largest over-threshold 
amplitudes and decreasing intermittency in the range below 
the first spectral kink. Several recent counter-intuitive obser- 
vational findings can be explained by the selective dissipation 
coupled to the nonlinear interaction among weakly dispersive 
KAWs. 



1 Introduction 

Kinetic Alfven waves (KAWs) are an extension of MHD 
Alfven waves in the range of high perpendicular wavenum- 
bers k± in the plane _L Bo, where linear and nonlinear effects 
due to finite values of p = k±pi become significant (Bo || .zis 
the background magnetic field, pi — Vn/^l% is the ion gyro- 
radius) (Hasegawa and Chen, 1976; Voitenko, 1998a). There 
are numerous observational and theoretical indications that 
MHD Alfven turbulence in the solar wind cascades towards 
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high k± and eventually reaches the KAW wavenumber range 
at the proton gyroradius scales, k±p. p ~ 1 (Leamon et al., 
1999; Bale et al., 2005; Alexandrova et al., 2008; Sahraoui 
et al., 2010). It is not yet certain what happens next with 
these KAWs: do they dissipate heating plasma (Leamon et 
al. 1999), or interact nonlinearly among themselves and pro- 
ceed cascading further towards higher k±, reaching electron 
scales (Alexandrova et al., 2008; Sahraoui et al., 2010). It 
has been envisaged that the nonlinear evolution and related 
wavenumber spectra in the range below the spectral break 
point k±i> ~ Pp 1 are dominated by MHD-type nonlinear in- 
teractions among Alfven waves, and the spectra at kj_ > k±_b 
are determined by linear and nonlinear properties of KAWs. 
If the dissipation is more efficient, the cascade should dissi- 
pate in the vicinity of k±i, and cannot reach electron scales, 
as was argued by Leamon et al. (1999), Howes et al. (2008), 
and Podesta et al. (2009) using Landau damping estimations. 
There are, however, observational indications (see Sahraoui 
et al., 2010, and references therein), that the nonlinear in- 
teraction among KAWs is faster than their dissipation, and 
the turbulence cascade proceeds further at higher k±p p 1 
creating a kinetic-scale turbulence of KAWs. 

From observational point of view, the transformation oc- 
curs at the spectral break points /& dividing shallower MHD 
spectra oc f~ 5 ^ 3 at / < /& and steeper kinetic spectra with 
power indexes ranging from —2 to —4 at / < ft, which 
are observed in the solar wind (/is the frequency in the 
satellite frame). Because of the large solar wind velocity, 
Vsw 3* Va, the Alfvenic time variations uja ~ k z VA are usu- 
ally much slower than the Doppler frequencies in satellite 
frame uid = k • Vsw (Tailor hypothesis). Then the satellite- 
frame frequency spectra are dominated by the Doppler fre- 
quency, 2-7t/ = \k z VA — k- Vsw\ ~|k-Vsty|, represent- 
ing wave-number spectra. As the solar wind turbulence is 
dominated by large perpendicular wave vectors k± 3> k z 
(Sahraoui et al., 2010; Luo and Wu, 2010, and references 
therein), satellites measure perpendicular wavenumber spec- 
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Fig. 1. Turbulent Alfvenic spectrum extending over three sub- 
ranges, with the steepest slope in the weakly/mildy dispersive KAW 
sub-range (interpretation of the spectrum shown in Fig. 1 by Chen 
et al., 2010) 



tra, / oc k±. There are of cause rare cases B || Vgy^, where 
frequency measures parallel wavenumber, / cx k z . The spec- 
tral break fi, of the turbulence is often associated with one of 
the proton kinetic scales, proton gyroradius p p or proton in- 
ertial length 5 p — Va/^u such that the observed break-point 
frequency 27r/& ~ Vsw / Pp or Vsw I $p (Leamon et al., 1999; 
Bale et al., 2005; Alexandrova et al., 2008; Sahraoui et al., 
2010). 

Because of the complex interplay between linear and non- 
linear dynamics of KAWs, theoretical interpretation of the 
turbulence, its dissipation, and related spectra in the KAW 
range is still incomplete. In particular, recent theoretical 
analysis by Podesta et al. (2010) argues that the KAW cas- 
cade subject to collisionless Landau damping cannot reach 
electron scales in the solar wind conditions, which contra- 
dicts the opposite conclusion by Sahraoui et al. (2010) based 
on observations. Using Cluster data, Sahraoui et al. (2010) 
found that the wavenumber spectra of MHD and KAW tur- 
bulences have slopes cx fcj 1 ' 7 and cx /cj 2 ' 8 , respectively, and 
the KAW turbulence extends to electron scales in the solar 
wind at 1 AU (Sahraoui et al., 2010). Between these fcj 1 7 
and fcj 2,8 spectra, Sahraoui et al. (2010) also noticed much 
steeper cx k J 4 spectra that appear in the weakly/mildy disper- 
sive KAW sub-range 0.6 < k±p p < 2 (see their Fig. 1). The 
same spectral form in the MHD/kinetic transition range, con- 
taining two spectral kinks with steepest spectra inbetween, 
can be seen in other recent studies - see example in Fig.l 
adopted from the paper by Chen et al. (2010) (see also Fig. 
1 by Smith et al., 2006). 

Steep variable spectra in the same wavenumber range were 
observed before (Leamon et al. 1999), but without connec- 
tion to shallower higher-wavenumber (higher-frequency in 



the satellite frame) spectra that were unavailable. These steep 
spectra were called the "dissipation range" spectra, and were 
associated with dissipation, mainly via kinetic ion-cyclotron 
and Landau damping. However, the nature of the "dissi- 
pation range" and its spectra is not so clear. For example, 
recent observations of reduced magnetic helicity imply the 
presence of counter-streaming KAWs surviving the "dissipa- 
tion range" rather than ion-cyclotron damping in it (Carbone 
et al. 2010). 

Analyzing ACE spacecraft data, Smith et al. (2006) have 
found that the larger spectral fluxes (as measured at 0.01 Hz) 
are followed by the steeper spectra in the "dissipation range" 
above the spectral break ~ 0.3 Hz. This counterintuitive 
observational fact is difficult to explain by ion-cyclotron and 
Landau damping. Smith et al. (2006) did not find any regu- 
lar dependence of on the cascade rate. However, later on 
Markovskii, Vasquez, and Smith (2008) studied the statistics 
and scaling of the spectral breaks and concluded that their 
positions are determined by a combination of their scales 
and the turbulent amplitudes at that scales, which suggests 
a non-linear dissipation mechanism for the solar wind tur- 
bulence. Again, kinetic ion-cyclotron and Landau damping 
mechanisms would not lead to such behavior of the dissipa- 
tion range. 

Motivated by these findings, in the present paper we ana- 
lyze the transition wavenumber range from MHD to KAW 
turbulence focusing on nonlinear KAW properties. We 
demonstrate that the observed spectral forms and steep spec- 
tra in the "dissipation range" can be explained by the nonlin- 
ear interaction of weakly dispersive KAWs without involving 
kinetic ion-cyclotron and Landau dissipation mechanisms. 



2 Weakly dispersive k±p p < 1 sub-range of the KAW 
turbulence 

In the weakly dispersive range, in the limit k\p^ <C 1, the 
rate of the nonlinear interaction among highly oblique k z <C 
k± co-propagating KAWs is calculated by maximizing the 
matrix element of the 3-wave KAW interaction calculated by 
Voitenko (1998a,b): 

Dl. 

(1) 

Tp -DO 

where B^ is the KAW amplitude at the (anisotropic) length 
scales X z — 2Tr/k z ; X± = 2n/k±, where k z and k± are 
parallel and perpendicular KAW wavenumbers. We put 
T e || /T p ± = 1 in (Q]) and further on for simplicity. 

In the KAW case k± k z and the nonlinear wave dynam- 
ics is driven by the energy exchange among short cross-field 
length scales, whereas the parallel scales follow the perpen- 
dicular ones kinematically (in the weak turbulence) or via 
critical balance (in the strong turbulence). Then the ampli- 
tude Bk can be related to the omnidirectional spectral energy 
density WkX as Bk = \fk~l_ Wk± (WLl is defined such that 
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J Q dk±Wk± = total fluctuation energy per unit volume). 
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The nonlinear interaction rate of counter-propagating /i <C 
1 KAWs is (Voitenko, 1998a,b) 



7fetl : 



:0.3Oi— n — . 

VTi -DO 



(2) 



In the weakly dispersive sub-range 7^ < 7^, but short 
linear correlation time among counter-propagating KAWs 
can reduce the counter-propagating nonlinear rate below 
7^. Note that there is no explicit k z - dependence in the 
above expressions for interaction rates. 

2.1 Weak KAW turbulence 

The nonlinear interaction among co-propagating KAWs can 
be considered weak if their nonlinear rate (H) is less than the 
dispersive part of the KAW frequency: < SuJk, where 
Soj k — k z VAH 2 in isothermal plasmas. In this case, the 
conservation law for the generalized enstrophy (dispersive 
part of energy) apply, and nonlinear interaction among co- 
propagating KAWs (0} establishes perpendicular wavenum- 
ber spectra 



W k oc k ± - 
W k oc kl 4 



(3) 
(4) 



created by the direct enstrophy and inverse energy cascades, 
respectively (Voitenko, 1998b). 

For the axially symmetric turbulence in the cross-field 
plane we can define a reduced omnidirectional spectral 
power W k ± = 2irk 1 _ J dk z B\, such that W = J dk±_W k ±. 
The energy exchange among different k± does not depend 
on k z (cf. eqs (6.1) and (6.2) by Voitenko (1998a)). Hence 
the weakly turbulent 3D spectra (|3]|4|i have their correspond- 
ing reduced omnidirectional power spectra 



W k ± oc k] 
W k ± oc k] 



(5) 
(6) 



The omnidirectional wavenumber spectra are those measured 
in the solar wind by satellites as ID Doppler frequency spec- 
tra if the solar wind velocity V$w -It Bo. In the case of tur- 
bulent spectra that are axially asymmetric around Bq axis, 
the measured ID spectrum may have a larger power index, 
approaching W k j_ oc fej in the extreme case of flat turbu- 
lence (say, oc fc,T 5 ). Accounting for a possible anisotropy, 
the steepest spectra produced by the weakly dispersive tur- 
bulence are 



Wk± ocfc_ L 4 ^fc J 



(7) 



Nonlinear interaction among counter-propagating KAWs 
(fJJ produces an imbalanced turbulence with weakly turbulent 
omnidirectional spectra 



Wkx^kJ 1 , (9) 

created by the direct enstrophy and inverse energy cascades, 
respectively (Voitenko, 1998b). The counter-propagating in- 
teractions cannot produce steep observed spectra in the tran- 
sition range. 

At first sight, interaction among counter-propagating 
KAWs in the weakly dispersive range is stronger than among 
co-propagating KAWs: 7^ 3> 7^ for However, 
short correlation times among counter-propagating KAWs, 
i~cx\, ~ A z /Va, can reduce their interaction strength as com- 
pared to the co-propagating KAWs that keep in phase longer 
time, r c ff ~ /j,~ 2 \ z /Va t c ^ for \l -C 1. In this case, the 
dominant (outward-propagating in the solar wind) compo- 
nent of the imbalanced turbulence will be shaped by the non- 
linear interaction among co-propagating KAWs with omni- 
directional spectra (0, as described above. 

The k z spectra of the weak KAW turbulence are deter- 
mined by the kinematics of three-wave resonant interactions: 

W k ^k- X ' 2 . (10) 

However, some non-kinematic factors, like finite resonance 
width, can make the parallel spectrum significantly different 
from ( [Tol l. This point needs further investigation. 

2.2 Strong KAW turbulence 

In the strongly turbulent regime the nonlinear time scale 
(~ l/7fe' i ) becomes equal to or larger than the linear 
one, (~ 1/ 5u>k), so that the enstrophy should not be con- 
served any more. In this regime we find the energy k±- 
spectrum from the condition that the energy flow through any 
k± does not depend on k±: 

e~ Bl/r^ L = const. (11) 

The nonlinear evolution time for co-propagating KAWs can 
be estimated as ~ 1/7^, where 7^ is given by ([T). 

Then from (fTTT i we find the scaling for the fluctuating mag- 
netic amplitude B k oc fcj 1 , which results in the omnidirec- 
tional energy spectrum 



W kA 



R2 

— ^oc£-~ 3 

fc, • 



(12) 



W k± ^kl 2 ; 



(8) 



Again, one can observe steeper spectra oc (fc^ 3 -f-fcj 4 ) if 
the strong KAW turbulence is not exactly axially symmetric 
around Bq. 

Since 7^ depends on k z only through B k , the k z - 
dependence can appear via any functional form with argu- 
ment involving some combination of k z and k±. Additional 
assumptions linking k z and k±, like critical balance hypoth- 
esis, will be studied in another paper. 

The strongly turbulent spectra of weakly dispersive 
counter-propagating KAWs can be found from ( fTTT i with 
l/rf i ~7^givenby©: 

B 2 

W kl _ ~ -± oc kl 7/3 (13) 
k± 
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3 Strongly dispersive k±p p > 1 sub-range of the KAW 
turbulence 

In the strongly dispersive sub-range of KAWs fj, > 1, 
which was considered in the literature and named as the 
"KAW range" (see Schekochihin et al., 2009, and refer- 
ences therein), the rate of nonlinear interaction among co- 
propagating KAWs is (Voitenko, 1998a;b) 



(14) 



VTi -DO 

For counter-propagating KAWs, the nonlinear interaction 
rate is almost the same, 



: 0.2^-^V|^. 

VTi -DO 



3.1 Weak turbulence < u k ) 



The weakly turbulent perpendicular wavenumber spectra of 
co-propagating KAWs, 



B 2 ocfc; 7/2 ; 
^cxfcj 3 , 



(15) 



are created by the direct energy and inverse enstrophy cas- 
cades, respectively (Voitenko, 1998b). Again, the nonlin- 
ear rate for the counter-propagating KAWs can be reduced 
by shorter linear correlation times as compared to the co- 
propagating KAWs. Therefore, the omnidirectional spectra 



W / fc _ L ocfc_ L 5/2 , 
W k ± oc fc7 2 



(16) 
(17) 



can be formed by strongly dispersive KAWs in the weakly 

— 5/2 

turbulent regime. Among these, the oc k , spectrum 
formed by the direct energy cascade is preferable. With lo- 
cal deviations from axially symmetry, one can expect steeper 
spectra oc /cj 2 ' 5 fcj/ 3 . 

3.2 Strong turbulence (uj k ~ j k L ) 

In the strong turbulence of co-propagating KAWs, the scaling 
of magnetic field amplitude B k with k± is found from the 
condition (Qj} where r k L ~ 1 /y k ^ with 7^ given by ( Q : 



Bk oc k 



-2/3 



This results in the familiar omnidirectional energy spectrum 
in k±: 



The "parallel" k z \\ Bq spectrum 

W k \\ oc fcj 2 

follows from the critical balance condition. 



(18) 



(19) 



4 MHD/kinetic Alfven transition 
4.1 Spectral kinks 

In the Goldreich and Sridhar (1995) MHD model, the AW 
nonlinear interaction rate at scale X± in the plane _L Bq can 
be written as 



7^ 



A 1 



-k ± V A ^. 
2n S ' 



(20) 



where v\± is the velocity and B k is the magnetic field am- 
plitude at the scale X± — 2ir/k^. The corresponding MHD 
AW spectrum B k oc kj 2 ^ 3 follows from the independence 
of the energy flux through k. The omnidirectional spectrum 
W k ~ B k /k±cc k ± °^ 3 is seen by satellites in the MHD range 
as the ID Doppler frequency spectrum. 

As the MHD and weakly dispersive KAW sub-ranges have 
very different slopes, the first spectral kink should appear 
at the wavenumber where their respective nonlinear interac- 
tion rates are equal. Comparing the nonlinear rates, 7^ = 
Tfcti' we obtain the spectral kink wavenumber /ii ~ 0.5, at 
which the ID spectrum should change from —(3/2-9-5/3) 
to — (3-j-4). The transition wavenumber for the 7j? s = 7^ 
transition is practically the same, /ii ~ 0.6. 

However, above estimations did not take into account 
the weakening of MHD nonlinear interactions by the dy- 
namic alignment between velocity and magnetic perturba- 
tions (Boldyrev, 2005) and/or by the nonlocal decorrelation 
mechanism proposed by Gogoberidze (2007). In general, the 
interaction rate can be written as a reduced GS rate (RGS) 



-.RGS p „,GS 
7fc - -«A_L7fc 



(21) 



with the scale-dependent reducing coefficient R k ±. Both 
Boldyrev's and Gogoberidze's phenomenologies give the 
same scaling for R k ±, 



R. 



fc_L ; 



VX± or A 1/4 



but with different normalization velocities Vjv, such that the 
Boldyrev/Gogoberidze ratio =vl/Va, where vl is the ve- 
locity amplitude at the driving scale L (wavenumber k^). 
Having in mind that the dynamic alignment saturates when 
approaching small scales, the actual value of the Gogob- 
eridze coefficient can be larger even in the case vl <Va- 
The reduced interaction rate proposed by Gogoberidze can 
be written as 



( 



k L 



-1/4 



2vr B 



(22) 



Given the typical width of the MHD inertial range in the so- 
lar wind kb±/ki J ~ 10 3 , we find that the interaction rate is 
reduced considerably in the vicinity of breakpoints, yj} GS ~ 
0.25 7fe GS 

As the nonlocal decorrelation mechanism implies counter- 
propagating MHD waves, the counter-propagating KAWs 
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should undergo the same decorrelation. But co-propagating 
KAWs do not suffer from such decorrelation, and therefore 
we consider here the MHD/kinetic transition dominated by 
the nonlinear interaction among co-propagating KAWs. In 
addition, the co-propagating KAWs can keep in phase much 
longer than the counter-propagating KAWs. We therefore 
use (Q~|i for kinetic and (fJTJ for MHD interaction rate, and es- 
timate the first spectral kink between the shallow MHD spec- 
tra — (3/2-9-5/3) and steep weakly dispersive KAW spectra 
-(3-5): 

^~0.6^/lhI. (23) 

With Gogoberidze's rate (|22| | /ii ~ 0.2. But one should bear 
in mind that there are a number of factors, including a par- 
tial turbulence dissipation, which contribute to Rki_ and can 
make it smaller or larger than the Gogoberidze's value. 

The second kink point should appear between weekly 
(Hp <C 1) and strongly (p p ^> 1) dispersive regimes of the 
KAW turbulence at 

M2>1, (24) 

where we allow for a possible building-up of still steeper 
slope just above /i = 1 in the cases where the MHD/KAW 
transition is not yet completed at k± = p~ l . The spectrum 
slope above /12, is —(2.5-5-3), which is significantly shal- 
lower than in the weakly dispersive range. 

4.2 Spectral forms 

The steepness of spectra in the weakly dispersive range de- 
pends on what kind of KAW turbulence picks up the turbu- 
lent cascade at /j, ~ 0.2, weak or strong. If the critical bal- 
ance condition holds at /x ~ /ii, then the turbulence of weakly 
dispersive KAWs is strong above In this case, strong 
KAW turbulence develops a steep energy spectrum oc fcj in 
the weakly dispersive sub-range, connecting shallower MHD 
(oc k ± °^ 3 ) and strongly dispersive KAW (oc k ± 7 ^ 3 ) spectra. 

Significantly steeper spectra in both KAW sub-ranges can 
be produced by the weak KAW turbulence and by local devi- 
ations from the azimuthal symmetry of the turbulence (up to 
about oc fcj] 4 ' 5 in the weakly dispersive sub-range, and oc kj 3 
in the strongly dispersive sub-range). Transition to the weak 
turbulence regime may be facilitated by the partial wave dis- 
sipation via non-adiabatic ion acceleration/stochastic heating 
that does not depend on k z but does depend on k± reducing 
larger amplitudes at k± > fcj_thr- In such a way, the critical 
balance between linear and nonlinear time scales is violated 
in favor of weak turbulent regime. After that, the weak turbu- 
lent cascade of KAWs develops above fcj_thr and establishes 
steepest KAW spectra. 

In both weak and strong turbulence regimes, the resulting 
spectra have two kinks, down and up, with the steepest slopes 
in between them in the weakly/mildy dispersive sub-range. 
In general, the "non-dissipative" scenario is as follows: the 




Frequency 



Fig. 2. Double-kink pattern produced by the MHD/weakly disper- 
sive/strongly dispersive turbulence transitions. The kinetic KAW 
spectra are given for the weakly turbulent regime (in the strongly 
turbulent regime, spectra -4 and -2.5 are replaced by -3 and -7/3, 
respectively) 



turbulence, driven at a large MHD scale (fedr — 27r/L(jr), 
develop the shallowest oc kj 3 ^ 2 kj 5 ^ 3 spectra in the MHD 
sub-range k^r <k±<k±i, then it proceeds as a KAW tur- 
bulence with steepest oc fcj 3 -j- /cj 4 ' 5 spectra in the weakly 
dispersive sub-range k±i <k±< k±2, and then above k±2 it 
proceeds as the KAW turbulence with oc k ± -j- k ± spec- 
trum in the strongly dispersive sub-range k±2 < k± < fcj_diss- 
This last sub-range may extend to the dissipative wavenum- 
ber fcj_diss at electron length scale (Sahraoui et al., 2010). 

Consequently, Alfvenic turbulent spectrum in the 
MHD/kinetic transition range attains an universal double- 
kink form (2-k pattern). This 2-k spectral pattern is shown 
schematically in Fig. 2 for the case of purely nonlinear 
non-dissipative transition. A (variable) slope of the weakly 
dispersive KAW spectrum depends on the ratio of turbulent 
energies cascading in strong and weak turbulent regimes, 
which can differ from case to case. The local slope should 
in principle lie between -3 and -4. But the shallower £-3 
spectra can be produced by a fraction of the MHD cascade 
extending above the first kink, and the steeper j-4 spectra 
can be produced by azimuthal asymmetry of the turbulence 
in the cross-field plane. 

A similar 2-k pattern can in principle be produced by the 
dissipative transition, suggested by many previous authors, 
but conditions required for that are rather special. Namely, 
the relative dissipation rate (as compared to the nonlinear in- 
teraction rate) should be much stronger in the range fcj^i < 
k± < k±2 than in the range k± > k±2- 

In any case, the presence of a high-wavenumber cascade 
and turbulence at k± > k± 2 imply the nonlinear transfer and 
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spectral flux in the range kj_i < k± < fcj_2 as well, which 
means the nonlinear interaction should be taken into account. 
In real situations the relative importance of the effects due 
to dissipation versus weak turbulence versus strong turbu- 
lence in the MHD/kinetic turbulence transition can be differ- 
ent from case to case. The 2-k pattern described above can 
be noticed in many high-resolution high-frequency Cluster 
measurements (see examples in papers by Kyiani et al., 2009; 
Sahraoui et al., 2010; Chen et al., 2010), and can be also no- 
ticed in some previous measurements where frequency ex- 
tended to 1 Hz or a little above (see e.g. Fig. 1 by Smith et 
al. 2006, showing ACE data). 

An actual wavenumber range where super-adiabatic ion 
acceleration and related wave damping come into play is also 
variable. The threshold behavior suggests that with stronger 
spectral fluxes it comes into play earlier and weakens the 
MHD turbulence. In the cases where the flatness follows the 
trends shown in Fig. 3 by Alexandrova et al. (2008) for Clus- 
ter data, the super-adiabatic ion acceleration and partial wave 
damping may be active well below the apparent spectral kink. 



5 Dissipation of KAWs 

In this section we introduce several pros and cons concerning 
basic dissipations mechanisms for KAWs, but their detailed 
investigation is postponed for future. 

Wu and Yang (2007) considered self-consistent velocities 
of minor ion species in KAW solitons and found them dis- 
tributed proportionally to the ion mass-to-charge ratio. How- 
ever, these velocities cannot be interpreted as thermal ones 
increasing temperature because they contribute to the non- 
thermal line broadening rather than thermal line width. A 
non-adiabatic disconnection from the wave fields is needed 
for the ions to gain some energy increase after the wave is 
passed by. Such process was considered by Voitenko and 
Goossens (2004), who shown that KAWs undergo strong 
non-adiabatic interaction with ions. This interaction require 
a certain threshold-like amplitude/wavelength relation for 
the dissipation switch-on, but does not require a long-time 
stochastic walk for the ions to gain a significant energy in- 
crease. 

Chandran et al. (2010) have shown that another process 
related to non-adiabaticity - stochastic plasma heating - can 
absorb up to half of the turbulent cascade flux at k±p p ~ 1. 
This result implies the MHD nonlinear rate, which may be 
not true at k±Pp ~ 1 where the KAW nonlinear interaction is 
faster and can pass more energy in the high-fcj^ range. 

Yet another nonlinear interaction of the broadband 
Alfvenic turbulence with ions, via nonlinear Landau damp- 
ing, was studied by Nariyuki et al. (2010), who shown that 
the ion heating proceeds both along and across the back- 
ground magnetic field and produces asymmetric ion velocity 
distributions. On the other hand, because of the quasi-linear 
platea formation in velocity distribution functions, classic 



Landau damping can be highly reduced in the weakly col- 
lisional solar wind (Voitenko and Goossens, 2006; Rudakov 
et al.,2011). 

5.1 Landau damping 

Parallel components of the KAW electric E z k and magnetic 
_B z k fields make the KAWs efficient in Cherenkov interac- 
tion with plasma species via kinetic mechanisms of Landau 
and transit-time damping that were commonly used in esti- 
mations of KAW dissipation. However, these mechanisms 
are based on the resonant wave-particle interactions that de- 
pend strongly on the local parallel slopes of particle veloc- 
ity distributions F s (V z ) at parallel velocities V z equal to the 
wave phase velocity Uk/k z . In particular, quasi-linear diffu- 
sion reduces resonant slopes and Landau damping (Voitenko 
2006): 



1L- 



1- 



TC 



TKAW 



(25) 



where jff is the Maxwellian Landau damping, and tkaw 
and tc are the characteristic diffusion times of parti- 
cles due to KAWs and Coulomb collisions, respectively. 
KAWs flatten F S (V Z ), Coulomb collisions restore it back 
to Maxwellian, and balance between two results in (l25T >. 
Whereas the detailed analysis of (l25T l as function of k± is 
quite complex (subject to separate study), our estimates, sim- 
ilar to that by Voitenko and Goossens (2006), show that for 
typical fluctuation level Wf ~ 10 _1 nT 2 /Hz at k±p p ~ 1 in 
the solar wind conditions tq /tkaw ^ 1 for both electrons 
and protons, and Landau damping is thus highly reduced. 
Therefore, conclusion by Podesta et al (2009) that the KAW 
turbulence cannot reach electron scales in the solar wind, 
based on the Maxwellian Landau damping, should by recon- 
sidered. 

5.2 Non-adiabatic threshold for turbulent dissipation 

The rate of the non-adiabatic cross-field acceleration of the 
ions i by oblique Alfven waves is (Voitenko and Goossens, 
2004): 

nAv A B ± v A )dxB 



— ibA 



1 



(26) 



where V iz is the parallel ion velocity, E± and B± are the 
Alfvenic electric and magnetic fluctuation, _L B^. 

Using E±/B± ~ Va/c in the weakly dispersive range, 
and ignoring possible field-aligned streaming of ions, the 
threshold-like condition for this kind of wave-particle inter- 
action, 



k±d p —> u 



(27) 



where z/,- = 0*/ Q, p is the threshold value for rjk above which 
the particular ion species i is heated super-adiabatically. 
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The condition for ion acceleration apply for any partic- 
ular ion specie, but the related wave dissipation depend on 
all ion species and their parameters, like species abundances, 
temperatures, etc. Nevertheless, condition for efficient wave 
dissipation can still be written in the form ( l27| > with a super- 
adiabaticity % in the left hand side, but with different thresh- 
old v w in the right hand side, which is not easy to find. One 
can guess that the wave threshold should be close to the ac- 
celeration threshold for the dominant ion species v w ~ z/j. 
Anyway, even without knowing the exact threshold value 
v w , it is possible to derive several useful scalings that can be 
tested observationally. So, for a power law scaling of mag- 
netic amplitudes, B\ oc k~J_ q , we obtain the spectral depen- 
dence of the super-adiabaticity rjk'- 

( fc, \ 1 ~ q/2 

Vk=mi \kh) ' (28) 

where Bki is the reference magnetic amplitude at the ref- 
erence length scale A_u = 2ir/kxi, and r)ki is the super- 
adiabaticity at k± = fcj_i- For the sake of convenience we 
choose the reference wavenumber equal to wavenumber of 
the first spectral kink k±i. 

Since q ~ 2/3 in the MHD range, super-adiabaticity rjk oc 

2 /3 

k"^_ grows with k± as long as k± < kxi- But the situa- 
tion is reversed in the weakly dispersive KAW range k± > 
k±i, where g~ 3 and super-adiabaticity decreases with k± 

— 1/2 

as rjk oc k ± . Such spectral k± -dependence of r)k indi- 
cates that the most favorable conditions for super-adiabatic 
ion heating and related wave dissipation are achieved in the 
vicinity of the first spectral kink, k± ~k±\. This is shown 
schematically in Fig. 3, where we used the omnidirectional 
spectral representation Wk oc B%/k± oc k^ with p = q+l. 
Since p < 2 (~ 5/3) in the MHD range, and p > 3 (~ 4) in 
the weakly dispersive KAW sub-range, the "threshold" spec- 
trum 

WWoc^cxfcX 3 , (29) 
k± 

that follows from the super-adiabatic condition, can fall be- 
low the observed turbulent spectrum Wk around the first 
spectral kink. 

Once the threshold rjk — v% is overcomed in some 
wavenumber range for some ion species, the ions enter 
regime of strong acceleration. In turn, because of its 
threshold-like character, the super-adiabatic ion heating pro- 
vides a highly selective dissipation mechanism for waves, 
affecting only strongest fluctuations with over-threshold am- 
plitudes. In principle, the ability of turbulence to produce 
intermittent large-amplitude fluctuations increases the value 
of driven parameter (l28t . where one should use the spectrum 
and amplitudes of intermittent fluctuations instead of the reg- 
ular turbulent spectrum. The eventual rate of the plasma 
heating and turbulence dissipation should follow from the 
balance between two processes: (i) production of the over- 
threshold intermittent fluctuations by the turbulence, and (ii) 




k±i k± 



Fig. 3. Typical Alfvenic turbulent spectrum (solid line) in the 
weakly/mildy dispersive KAW sub-range and "threshold" turbulent 
spectrum (dashed line) required for the super-adiabatic ion accel- 
eration. Super-adiabatic acceleration is possible around the first 
spectral kink, where the turbulent spectral power raises above the 
threshold spectral power. 



accommodation of turbulence energy by accelerated ions and 
its further redistribution into the bulk plasma. Helios ob- 
servations have shown that the flatness (measure of inter- 
mittency) increases with increasing wavenumber (Alexan- 
drova et al., 2008), which progressively increases the super- 
adiabaticity parameter above the value given by (l28t . Then, 
at some large enough wavenumber, the level of intermittent 
fluctuations can reach the super-adiabatic threshold even if 
the regular turbulent level remain below it. . 

Dissipation due to super-adiabatic heating/acceleration 
tends to reduce the over-threshold fluctuations at every scale 
to the corresponding threshold value given by (1271 . Then, 
in accordance to ( |27] |. the upper bound for the reduced in- 
termittent amplitudes scales as J3 2 hr oc fcj_ 2 , and since the 
magnetic power spectrum in this range has shallower scal- 
ing B\ oc fc ± 2 ^ 3 , the flatness (and higher order normalized 
structure functions as well) should decrease with wavenum- 
ber in the MHD range below the first spectral kink. This 
can explain another interesting feature, a local decrease of 
the flatness in the spacecraft frequency range 0.02-^0.2 Hz 
(which is still below the apparent spectral kink), found by 
Alexandrova et al. (2008) using Cluster data. We suggest 
that such behavior of the flatness may indicate a partial dissi- 
pation of Alfven waves via super-adiabatic ion acceleration 
in the corresponding wavenumber range. 

In turn, highly anisotropic ion distributions are produced 
by super-adiabatic acceleration (Voitenko and Goossens, 
2004), which can drive anisotropic ion-cyclotron instabili- 
ties redistributing energy further. Since the super-adiabatic 
ion acceleration is very fast, within a fraction of the cor- 
responding ion gyroperiod, the quasi-stationary rate of tur- 
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bulent dissipation will be determined by the ion-cyclotron 
instability increment. The situation is thus more complex 
here and opposite to that observed in hydrodynamics, where 
viscosity washes out smallest amplitudes when approaching 
the dissipation range, but large-amplitude fluctuations sur- 
vive increasing intermittency. The behavior of the intermit- 
tency found by Alexandrova et al. (2008) is not typical for 
the linear Landau damping as well. 

After the relative perpendicular/parallel power in the spec- 
trum and the strength of MHD interaction are reduced, the 
transition to the weak KAW turbulence is made possible and 
leads to the steepest spectra in the weakly dispersive sub- 
range. It would be interesting to analyze intermittency by the 
Rank-Ordered Multifractal Analysis (ROMA) (see Chang et 
al., 2010, and references therein), which would allow to find 
out if the solar wind turbulence possesses different fractal 
properties in the three mentioned above sub-ranges. 

6 Summary and Discussion 

For the first time, a weakly dispersive sub-range of the 
KAW turbulence is distinguished and studied in the con- 
text of MHD/kinetic turbulence transition. We show that 
the KAW turbulence and its spectra in the weakly dispersive 
sub-range differ significantly not only from the conventional 
MHD Alfvenic turbulence, but also from the strongly dis- 
persive KAW turbulence. Namely, the nonlinear interaction 
of weakly dispersive KAWs is capable to produce steepest 
spectra oc fcj 3 /cj 5 in the wavenumber range fcxi < k±_ < 
p~ l , connecting shallow MHD spectra oc k ± 3 ^ 2 ~ k ± J ^ 3 be- 
low the first spectral kink, k± < fcj_i, and "intermediate" 

7/3 o 

oc k ± 1 -rfcjj spectra of strongly dispersive KAWs above 
the second spectral kink, k±> p^ 1 . 

The universal spectral form - 2-k pattern resulting from 
such spectral dynamics in the transition range is shown 
schematically in Fig. 2. Turbulent spectra observed recently 
by the Cluster spacecraft often exhibit such 2-k pattern in the 
transition wavenumber range (see for example Fig.l). It is 
still not certain what is the role of Landau damping in pro- 
ducing so steep spectral drop at k±p p < 1. Any kind of ki- 
netic dissipation in the weakly collisional solar wind should 
be self-consistently saturated at a reduced level by the lo- 
cal platea formation in the velocity distribution functions of 
plasma species. At least a quasi-linear theory is needed to ac- 
count for the particles' feedback reaction on the energy input 
from the waves, and numerous previous estimations based on 
the Maxwellian Landau damping should be re-evaluated. 

Podesta (2009) reported a significant flattening of the high- 
frequency parallel spectra and suggested it may be due to a 
plasma instability injecting a fraction of parallel propagating 
waves. On the other hand, this flattening can be produced by 
the transition to the weak KAW turbulence, possessing (in an 
ideal case) a very shallow spectrum (ITOb . However, because 
of many interfering factors, it is not certain if the parallel 



wavenumber spectrum (TTOb can be realized in the solar wind. 
The perpendicular wavenumber spectra are determined by 
the nonlinear interaction among perpendicular length scales 
and are thus quite robust. But the corresponding parallel 
wavenumber dynamics and spectra follow the perpendicular 
wavenumber dynamics and are often defined from a suitable 
functional form linking them to the perpendicular ones. This 
functional form may depend on a number of factors, includ- 
ing strength of the turbulence, partial turbulence dissipation, 
etc. In the extreme cases of weak and strong turbulence, the 
parallel dynamics is fixed, respectively, by the perpendicular 
one kinematically (via resonant conditions) and by adjusting 
linear and nonlinear time scales (via critical balance condi- 
tion). 

One can expect a high variability of spectral slopes in 
the weakly dispersive wavenumber sub-range, resulting from 
a mixture of several "pure" spectra that can be produced 
by KAWs in this range. In addition, our analysis suggests 
the super-adiabatic and/or stochastic cross-field acceleration 
of the solar wind ions as feasible mechanisms for a partial 
KAWs dissipation operating in the vicinity of first spectral 
kink. Both these mechanisms share the same non-adiabatic 
threshold and imply a selective dissipation of the over- 
threshold fluctuations with largest amplitudes. This kind of 
dissipation reduces high-amplitude intermittent fluctuations 
and should therefore produce local decrease of flatness in the 
dissipation range. Albeit there are observational indications 
for such behavior of flatness (see Fig. 3 by Alexandrova et 
al., 2008), this point needs further observational support. 

It seems that the synergetic action of selective wave dis- 
sipation and weak turbulence of KAWs influences both the 
spectral kink positions and the spectral slopes, making them 
dependent of the turbulence level. Namely, rjk, product of 
the turbulent amplitude and corresponding wavenumber, is 
the parameter facilitating transition to the weak KAW tur- 
bulence with its steeper spectra. As the spectral flux ~ rjf., 
the larger spectral fluxes imply larger rjk, which in turn im- 
ply steeper spectra in the weakly dissipative range. Such a 
counter-intuitive trend was found by Smith et al. (2006). 

On the other hand, in the vicinity of spectral kinks the non- 
adiabatic wave-particle interaction tends to reduce rjk to a 
near-threshold value, which results in the scaling B^i ~ fcj*. 
This scaling offer an explanation for observed spectral kink 
wavenumbers that were found to be inversely proportional 
to the fluctuation amplitudes at spectral kinks (Markovskii et 
al., 2008). 

Contrary to MHD Alfven waves, the dispersion law of 
KAWs, even weakly dispersive, is not degenerated with re- 
spect to k±. This makes possible 3-wave interactions with 
all 3 waves residing on the KAW branch, and there is no 
need in a zero-feii, k± ^ mode mediating the MHD turbu- 
lent cascade. Consequently, an additional spectrum of the 
KAW turbulence can be created by the cascading enstrophy 
(dispersive part of energy). The energy and the enstrophy 
cascade in opposite directions from the injection wavenum- 
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ber. As the turbulence of KAWs in the solar wind is driven 
at largest near-MHD length scales, it naturally proceed to 
smaller scales following a direct cascade route. In other cir- 
cumstances, and with different positions of the driving scale, 
one may observe inverse (e.g. Lui et al., 2008), or dual spec- 
tral transport, which is not easy to discriminate and describe 
in terms of cascades because of the non-local contaminations 
and scale mixing in a real finite-size and high-variable envi- 
ronment (see Voros et al., 2010). 

Again, contrary to the MHD AW turbulence, the 
KAW turbulence does not require the pre-existing counter- 
propagating waves for efficient cascading. Nonlinear inter- 
action among co -propagating KAWs is strong enough to es- 
tablish a co-propagating (completely imbalanced) KAW tur- 
bulence without involving the counter-propagating KAWs. 
If the co-propagating KAW turbulence develops in some 
wavenumber range (e.g. at k±i <k^< p" 1 ), then the ratio 
of sunward/anti-sunward Poynting fluxes should be frozen 
and remain approximately constant at these wavenumbers. 
This would provide another observational benchmark for the 
KAW turbulence, but we are not aware about such observa- 
tions so far. 
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